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1  Introduction 


Background 

Contamination  to  air,  surface  water,  and  groundwater  from  military- 
related  compounds  (MRCs)  has  been  of  increasing  concern  in  the  past  sev¬ 
eral  years  due  to  the  closing  of  munition  plants  and  bases.  MRCs  usually 
refer  to  some  type  of  explosive  compound  or  to  a  by-product  of  such. 
However,  some  radioactive  materials  may  also  be  classified  as  MRCs. 

Because  of  the  many  mediums,  or  multimedia,  available  for  the  spread 
of  contamination  (air,  surface  water,  groundwater,  etc.),  the  effects  of  the 
contamination  on  ecological  units  from  each  of  the  primary  modes  of 
exposure  need  to  be  evaluated  via  the  specific  protocol  of  the  exposure 
routes. 

The  assessment  of  the  exposure  to  ecological  units  from  the  many  medi¬ 
ums  through  the  use  of  an  exposure  model  is  desirable  to  determine  the 
risk  associated  with  exposure  to  MRCs  from  these  media.  This  assessment 
is  valuable  for  permitting  and  planning  activities  as  well  as  for  possible 
cleanup  operations  of  contaminated  sites  should  the  risk  to  an  ecological 
unit  become  too  great. 

Ideally,  when  a  general  idea  of  the  level  of  risk  associated  with  an 
exposure  is  desired,  the  use  of  a  screening-level  model  will  maximize  the 
amount  of  information  provided  while  minimizing  the  amount  of  effort 
required  to  obtain  the  necessary  information  to  make  an  assessment  of 
risk.  A  screening-level  model  refers  to  the  use  of  simplified,  quantitative, 
predictive  methods  that  minimize  time  and  effort  for  implementation.  Sim¬ 
plification  is  achieved  by  making  assumptions  that  reduce  complexity  of 
the  predictive  mathematical  formulations  and  input  data  requirements. 

The  intent  here  was  to  evaluate  the  existing  multimedia  exposure 
assessment  models  and  determine  which  models  would  lend  themselves 
easily  to  modification  to  handle  the  assessment  of  ecological  risk  of  expo¬ 
sure  from  military  unique  compounds.  The  primary  MRCs  addressed  in  this 
report  are  2,4,6-trinitrotoluene  (TNT),  hexahydro-l,3,5-trinitro-l,3,5- 
triazine  (RDX),  and  octahydro-l,3,5,7-tetranitro-l,3,5,7-tetrazocine 
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(HMX).  Also  in  this  report  is  a  brief  section  on  radioactive  contaminants 
such  as  plutonium,  americium,  etc. 


Objectives 

The  objectives  of  this  report  are  to  provide  a  review  of  the  screening- 
level  risk  assessment  models  evaluated  in  this  study  and  to  discuss  the 
modifications  of  the  selected  model  in  order  to  allow  it  to  address  the 
assessment  of  ecological  risk  of  exposure  to  military  unique  compounds. 


Chapter  1  Introduction 


2  Review  of  Exposure 
Assessment  Models 


A  review  of  existing  exposure  assessment  models  was  conducted  to 
meet  the  requirements  of  the  Fate  and  Effects  Research  Program  for  a 
screening-level  multimedia  exposure  assessment  model  for  MRCs.  A 
literature  research  effort  was  conducted  to  determine  the  most  useful 
exposure  assessment  models  for  this  task.  The  research  turned  up  numer¬ 
ous  models  which  were  narrowed  to  nine  potential  candidates  based  on 
their  use  as  a  screening-level  tool.  The  models  that  were  ultimately  consid¬ 
ered  as  potential  candidates  are  shown  in  Table  1. 


Table  1 

List  of  Exposure  Assessment  Models  Evaluated 

MEPAS  (Buck  etal.  1995) _ 

FRAMES1 _ 

MULTIMED  (Health  and  Welfare  Canada  1992) _ 

MULTIMED-DP  (Salhotra  et  al.  1 993) _ 

RESRAD-CHEM  (EPA  1995) _ _ _ 

MMSOILS  (EPA1995) _ _ _ 

PIRANHA  (Burns  1992) _ 

PCGEMS  (Access  EPA  1993) _ 

PRECIS  (Sandia  National  Laboratories  1994) 

1  Personal  Communication,  1997,  Gene  Whelan,  Pacific  Northwest  National  Laboratory,  Richland,  WA. _ 


Model  Histories  and  Backgrounds 


MEPAS 

The  Multimedia  Environmental  Pollutant  Assessment  System  (MEPAS) 
was  developed  by  Pacific  Northwest  Laboratory  for  the  U.S.  Department 
of  Energy  (DOE)  Office  of  Environment,  Safety,  and  Health.  It  was 
derived  from  the  Remedial  Action  Priority  System  (RAPS)  and  is  com¬ 
prised  of  a  set  of  computer  programs  used  to  prioritize  hazardous, 
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radioactive,  and  mixed-waste  sites  based  on  their  potential  impact  on 
human  health.  MEPAS  is  applicable  to  a  wide  range  of  environmental 
management  and  regulatory  applications.  The  model  can  be  applied  to 
inactive  as  well  as  to  active  air  and  water  release  sites.  The  model  itself 
uses  a  simplified  contaminant-transport  and  exposure  algorithm  encapsu¬ 
lated  into  a  user-friendly  shell  containing  an  extensive  database  of 
approximately  476  contaminants.  MEPAS  is  an  integrated  system  of  ana¬ 
lytical,  semi-  analytical,  and  empirically  based  mathematical  models  that 
simulate  the  release  of  contaminants  from  sources,  migration  and  fate  of 
these  contaminants  through  various  environmental  media,  concentrations 
at  designated  receptor  locations,  and  health  effects  to  specified  sensitive 
receptors  caused  by  exposure  to  these  contaminants.  It  was  designed  to 
address  long-term  average  environmental  conditions  resulting  from 
release  to  the  environment  from  point  and  area  sources.  MEPAS  does  not 
address  acute  exposures  to  environmental  contamination. 


FRAMES 

The  Framework  for  Risk  Analysis  and  Management  of  Environmental 
Systems  (FRAMES)  was  developed  by  Pacific  Northwest  National  Labora¬ 
tory  for  the  DOE.  FRAMES  is  an  object-oriented  model  that  is  still  under 
development.  Within  FRAMES  will  reside  a  collection  of  computer  algo¬ 
rithms  that  will  simulate  the  following  elements  of  a  transport,  exposure, 
and  risk-assessment  system:  contaminant  source  and  release  to  environ¬ 
ment  (including  surface  hydrology),  overland  transport,  vadose-zone 
transport,  saturated-zone  transport,  atmospheric  transport,  surface  water 
transport,  food-supply  transport  (including  animals  and  plants  to  humans), 
intake  computation,  and  health  impacts.  In  FRAMES,  each  of  the  above- 
listed  elements,  and  many  that  are  not  listed,  can  and  will  be  represented 
by  separate  modules.  FRAMES  utilizes  three  types  of  files:  primary  data 
communication  files  which  transfer  information  between  modules,  global 
input  data  files  where  all  user  input  is  stored,  and  other  data  files  which 
allow  for  imported  data  files,  exported  data  files,  and  maintained 
databases. 

Within  FRAMES,  data  transfer  linkages  are  designed  to  avoid  data 
redundancy,  to  provide  dynamic  linked  modules  (modules  are  linked  as 
the  direct  result  of  user  selection),  and  to  provide  a  framework/user  inter¬ 
face  that  allows  for  relatively  easy  inclusion  of  additional  modules  and 
models,  promotes  access  to  national  databases,  minimizes  data-exchange 
requirements,  and  allows  for  unlimited  access  to  data. 


MULTIMED 

The  Multimedia  Exposure  Assessment  Model  (MULTIMED)  was  devel¬ 
oped  by  the  Center  for  Exposure  Assessment  Modeling  (CEAM)  of  the 
U.S.  Environmental  Protection  Agency  (USEPA).  It  simulates  the  move¬ 
ment  of  contaminants  leaching  from  a  waste  disposal  facility  or  from 
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contaminated  soils.  It  is  comprised  of  several  modules  that  predict  concen¬ 
trations  at  a  receptor  due  to  transport  in  the  subsurface,  surface  water,  or 
air.  The  model  is  accompanied  by  a  preprocessor,  PREMED,  and  a  post¬ 
processor,  POSTMED.  The  preprocessor  guides  the  user  in  the  creation  of 
a  Subtitle  D  input  file,  while  the  postprocessor  is  used  to  generate  plots  of 
data  obtained  from  the  simulation. 

The  model  calculates  contaminant  concentrations  in  groundwater,  sur¬ 
face  water,  and  air  resulting  from  emissions  from  hazardous  waste  sites. 
The  hydrologic  and  atmospheric  pathways  are  addressed  independently 
with  no  transfer  of  contaminant  between  the  two  pathways. 


MULTIMED-DP 

The  Multimedia  Exposure  Assessment  Model  -  Daughter  Products 
(MULTIMED-DP)  is  essentially  the  same  as  the  MULTIMED  except  that 
it  can  be  used  to  simulate  the  transport  and  fate  of  first-  and  second- 
generation  transformation  (daughter)  products  that  migrate  from  a  waste 
source  through  the  unsaturated  and  saturated  zones  to  a  down  gradient 
receptor  well. 


RESRAD-CHEM 

The  Residual  Radioactivity  (RESRAD)  model  was  developed  in  the 
early  1980s  by  Argonne  National  Laboratory  for  the  DOE.  It  was  initially 
developed  to  calculate  site-specific  residual  radioactive  material  guide¬ 
lines  and  radiation  dose/risk  to  an  onsite  individual  at  a  radioactively  con¬ 
taminated  site.  The  RESRAD  code  has  been  continuously  updated  and 
improved  to  now  include  chemical  contaminants  (RESRAD-CHEM). 

The  pathways  considered  by  this  model  include  direct  external  expo¬ 
sure,  dust  inhalation,  exposure  to  radon,  and  ingestion  of  plant  foods, 
meat,  milk,  aquatic  foods,  water  and  soil.  The  model  has  a  built-in  capabil¬ 
ity  to  do  sensitivity  analysis  on  input  parameters,  and  it  also  has  the  capa¬ 
bility  to  do  uncertainty  analysis  on  input  parameters  and  databases  using 
the  Latin  Hypercube  Sampling  (LHS)  technique. 


MMSOILS 

The  MMSOILS  model  was  developed  by  Tetra  Tech,  Inc.,  for  the 
USEPA.  MMSOILS  is  a  methodology  for  estimating  the  human  exposure 
and  health  risk  associated  with  releases  of  contamination  from  hazardous 
waste  sites.  The  methodology  addresses  the  transport  of  a  chemical 
through  groundwater,  surface  water,  soil,  the  atmosphere,  and  aquatic/ 
terrestrial  food  chains.  For  multimedia  exposure,  the  model  provides  esti¬ 
mates  of  human  exposure  through  individual  pathways  and  cumulative 
exposure  through  all  pathways  considered.  The  risk  associated  with  the 
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total  exposure  dose  is  calculated  based  on  the  chemical-specific  toxicity 
data. 

The  human  exposure  pathways  considered  include  soil  ingestion,  air 
inhalation  of  volatile  compounds  and  particulate  matter,  dermal  contact, 
ingestion  of  drinking  water,  consumption  of  fish,  consumption  of  plants 
grown  on  contaminated  soil,  and  consumption  of  animals  that  grazed  on 
contaminated  pastures.  For  multimedia  exposure,  the  methodology  pro¬ 
vides  estimates  of  human  exposure  through  individual  pathways  and  cumu¬ 
lative  exposure  through  all  pathways  considered. 


PIRANHA 

The  Pesticide  and  Industrial  chemical  Risk  ANalysis  and  Hazard 
Assessment  (PIRANHA)  model  was  developed  by  the  USEPA.  PIRANHA 
is  a  computer-based  system  for  applying  the  ecotoxicological  and  environ¬ 
mental  sciences  to  ecological  risk  assessment.  PIRANHA  includes  a  Geo¬ 
graphic  Information  System  (GIS-)  based  facility  for  locating  biological 
resources  potentially  at  risk  from  pesticide  and  industrial  chemicals,  a 
toxicological  inferencing  program  for  species-to-species  extrapolation  of 
acute  toxicity,  models  for  predicting  ambient  chemical  exposures  in 
aquatic  and  agro-ecosystems,  and  bioaccumulation  models  for  fauna. 

The  PIRANHA  model  can  use  up  to  three  sub-models  during  the  course 
of  analysis.  These  sub-models  are  the  Pesticide  and  Root  Zone  Model 
(PRZM),  the  Exposure  Analysis  Modeling  System  (EXAMS),  and  the 
Food  and  Gill  Exchange  of  Toxic  Substances  (FGETS)  model. 

The  PRZM  is  a  one-dimensional,  dynamic,  compartmental  model  for 
simulating  vertical  chemical  movement  in  unsaturated  soil  systems  within 
and  immediately  below  the  plant  root  zone.  It  allows  the  user  to  simulate 
movement  of  potentially  toxic  chemicals,  particularly  pesticides,  that  are 
applied  to  the  soil  or  to  plant  foliage.  Dynamic  simulations  allow  for  the 
consideration  of  pulse  loads,  the  prediction  of  peak  events,  and  the  estima¬ 
tion  of  time-varying  mass  emission  or  concentration  profiles. 

The  EXAMS  model  is  an  interactive  computer  program  intended  to 
give  decision-makers  access  to  a  responsive,  general,  and  controllable 
tool  for  readily  deriving  and  evaluating  the  behavior  of  synthetic  chemi¬ 
cals  in  the  environment.  The  model  estimates  exposure,  fate,  and  persist¬ 
ence  following  release  of  an  organic  chemical  into  an  aquatic  ecosystem. 

The  FGETS  model  simulates  the  bioaccumulation  of  nonmetabolized 
organic  chemicals  in  fish.  The  model  is  formulated  in  such  a  way  that  its 
parameterization  does  not  rely  upon  specific  toxicokinetic  studies  but, 
rather,  depends  on  physical  and  chemical  properties  of  the  compounds 
that  can  be  estimated  by  other  models  and  on  ecological,  morphological, 
and  physiological  characteristics  of  the  fish  that  can  be  obtained  from 
computerized  databases. 
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PCGEMS 


The  Graphical  Exposure  Modeling  System  PC  (PCGEMS)  is  a  stand¬ 
alone  version  of  the  Graphical  Exposure  Modeling  System  (GEMS)  for 
use  on  a  PC.  GEMS  began  development  in  1981  by  the  USEPA.  PCGEMS 
supports  exposure  and  risk  assessments  by  providing  access  to  single 
medium  and  multimedia  fate  and  exposure  models,  physical-chemical 
property  estimation  techniques,  statistical  analysis,  graphics  and  mapping 
programs  with  related  data  on  environments,  sources,  receptors,  and  popu¬ 
lations.  PCGEMS  is  an  easily  learned  interface  to  various  models,  pro¬ 
grams,  and  data  needed  for  exposure  and  risk  assessments.  Currently, 
PCGEMS  contains  associated  environmental  and  1980  population  data 
encompassing  most  of  the  United  States. 

The  PCGEMS  allows  the  use  of  mathematical  algorithms  to  calculate 
an  estimate  of  pollutant  concentration  based  on  a  number  of  environ¬ 
mental  factors.  It  also  allows  the  use  of  several  screening-level  models 
that  predict  chemical  partitioning,  the  separation  of  components  of  a 
chemical  mixture,  among  the  environmental  media.  These  models  include 
media-specific  capabilities  for  modeling  pollutant  releases  to  the  atmos¬ 
phere,  surface  water,  soil,  and  groundwater. 


PRECIS 

The  Probabilistic  Risk  Evaluation  and  Characterization  Investigation 
System  (PRECIS)  was  developed  by  Sandia  National  Laboratory  for  the 
DOE.  The  software  aids  site  investigators  in  performing  probabilistic  risk 
assessments  to  quantify  risk  to  human  and  environmental  receptors  (with 
associated  uncertainty)  resulting  from  site  contamination,  in  determining 
risk-based  action  or  cleanup  levels,  and  in  prioritizing  data  needs  for  char¬ 
acterization.  It  can  be  used  to  estimate  dose  and  action  or  cleanup  levels 
for  radionuclides,  risks  and  toxicity  effects  from  hazardous  chemicals, 
and  action  or  cleanup  levels  for  hazardous  chemicals.  The  software  incor¬ 
porates  a  user-friendly  Graphical  User  Interface  (GUI)  and  is  available  in 
versions  for  Windows  3.1®,  Windows  NT®,  Windows  95®,  and 
Macintosh®. 


Exposure  Model  Pathways  and  Features 

The  candidate  models  have  various  features  and  pathways  which  are 
summarized  in  Table  2.  A  more  detailed  comparison  of  the  models,  which 
includes  comparison  of  such  features  as  uncertainty  analysis,  partitioning, 
ability  to  simulate  transformation  products,  etc.,  is  shown  in  Table  3. 
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Table  2 

Pathway/Feature  Comparison  of  Candidate  Multimedia  Models 

Pathway/ 

Feature 

Model 

PCGEMS 

MEPAS 

MULTIMED 

m 

Transport  Media 

mg 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

N 

Y 

Y 

Y 

Soil 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Surface  water 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

N 

Y 

Source  Type 

Underground 
storage  tank 

N 

Y 

N 

N 

Y 

N 

Y 

N 

i 

Surface  water 

Y 

Y 

N 

N 

Y 

N 

Y 

N 

i 

Y 

Y 

Y 

Y 

Y 

N 

Y 

Y 

1 

Stack 

Y 

Y 

N 

N 

N 

N 

Y 

N 

1 

Soil 

Y 

Y 

N 

N 

Y 

Y 

Y 

Y 

i 

Road 

N 

N 

N 

N 

N 

N 

N 

N 

i 

Genera!  Features 

User  interface 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Default  data 

Y 

Y 

N 

N 

Y 

Y 

Y 

Y 

Risk 

assessment 

Y 

Y 

N 

N 

Y 

Y 

Y 

Y 

Y 

Toxics 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Nutrients 

N 

N 

N 

N 

N 

N 

N 

N 

N 

GW  Pathways 

Model  type 

Flow/Trans 

Flow/Trans 

Flow/Trans 

Flow/Trans 

Flow/Trans 

N/A 

Flow/Trans 

Flow/T  rans 

Transport 

Zone 

Both 

Both 

Both 

Both 

Both 

N/A 

Both 

Both 

Both 

Degradation 

Y 

Y 

Y 

Y 

Y 

N/A 

Y 

Y 

_ 

Retardation 

Y 

Y 

Y 

Y 

Y 

N/A 

Y 

Y 

_ 

Dispersion 

Y 

Y 

Y 

Y 

Y 

N/A 

Y 

N 

_ 

NPS  Pathways 

Watershed 

Agricultural/ 

Undeveloped 

Agricultural/ 

Landfill 

Landfill 

Landfill 

Landfill 

Agricultural 

Agricultural/ 

Landfill 

Landfill 

Agricultural 

Watersheds 

Single 

Single 

Single 

Single 

!■ 

Events  . 

Continuous 

mM 

BNI 

Continuous 

Continuous 

— 

Annual 

average 

Annual 

average 

Sediment 

transport 

Y 

Y 

N 

N 

IIHI 

Y 

Y 

N 

Adsorption 

Y 

N 

N 

N 

N 

N 

N 

Y 

Volatilization 

Y 

Y 

Y 

Y 

N 

N 

N  | 

i _ 

Continued) 
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Table  2  (Concluded) 

Model 

Pathway/ 

Feature 

PCGEMS 

MEPAS 

MULTIMED 

MULTI  MED- 
DP 

MMSOILS 

PIRANHA 

FRAMES 

RESRAD- 

CHEM 

PRECIS 

SW  Pathways 

Water  body 
type 

Lake/River/ 

Estuary 

River 

River 

River 

Lake/  River 

Lake/River 

River 

N/A 

Stream/ 

Pond 

Degradation 

Y 

N 

Y 

Y 

Y 

Y 

N 

N/A 

— 

Adsorption 

Y 

N 

Y 

Y 

Y 

Y 

N 

N/A 

— 

Volatilization 

Y 

N 

Y 

Y 

Y 

Y 

N 

N/A 

N 

Air  Pathways 

Release 

mechanism 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

N 

Dilution/ 

transport 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

— 

Washout 

N 

Y 

N 

N 

N 

N 

Y 

N 

N 

Deposition 

Y 

Y 

N 

N 

Y 

Y 

Y 

Y 

— 

Complex  terrain 

Y 

Y 

N 

N 

N 

N 

Y 

N 

N 

Wake  effects 

Y 

Y 

N 

N 

N 

N 

Y 

Y 

N 

Degradation 

Y 

Y 

Y 

Y 

N 

N 

Y 

Y 

— 

Note:  GW  =  Gro 
1  Calculations  ar 

undwater,  NPS  =  Non-Point  Source,  SW  =  Surface  Water 
e  based  on  time  since  material  was  placed. 

Selection  of  Screening-Level  Model  for 
Simulating  MRCs 

Based  on  a  model’s  overall  capability  as  a  screening-level  tool  for 
evaluating  MRCs,  the  ease  with  which  adaptations  or  modules  could  be 
made  to  an  existing  model  without  significant  code  modification,  and  the 
relative  ease-of-use  or  user-friendliness  of  the  final  product,  risk  assess¬ 
ment  models  were  chosen  as  the  building  blocks  for  the  creation  of  a  mod¬ 
ule  which  contains  the  primary  process  descriptions  for  the  fate  of  target 
MRCs.  The  models  that  were  chosen  for  the  adaptation  to  simulate  MRCs 
were  FRAMES  and  MEPAS. 


Evaluation  of  MEPAS  database  for  existing  MRCs 

The  MEPAS  risk  assessment  model’s  database  was  analyzed  to  deter¬ 
mine  the  current  chemical  property/risk  association  values  of  any  existing 
MRCs.  The  MEPAS  database  contained  information  for  TNT,  RDX,  and 
HMX.  The  existing  information  contained  in  the  MEPAS  database  on  each 
of  these  three  chemicals  is  presented  in  Table  4. 
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Table  3 

Detailed  Comparison  of  Potentially  Useful  Multimedia  Models 

Model 

Property 

MEPAS 

MULTIMED 

MULTIMED- 

DP 

RESRAD- 

CHEM 

MMSOILS 

FRAMES 

PIRANHA 

Uncertainty 

estimation 

N  (Sum3  is 
available 
separately  as 
option) 

Y  (Monte  Carlo) 

Y  (Monte  Carlo) 

Y 

Y 

Y 

Y 

Screening  level 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Assessment  level 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Risk  estimation 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Applicable  to 
NAPLs1 

N 

N 

N 

N 

N 

N 

N 

Homogeneous 

medium 

Y 

Y 

Y 

Y 

Y 

Y 

N/A 

Heterogeneous 

medium 

N 

N 

N 

N 

N 

N 

N/A 

Overland 

transport 

Y  (modified 

curve-number 

technique) 

N 

N 

N 

Y 

Y 

Y 

Sediment 

transport 

Y  (Universal  Soil 
Loss  Eqn.) 

N 

N 

Y 

Y 

Y 

N 

Evapo- 

transpiration 

Y  (Penman 
Method  with 
Correction 

Factor) 

Y 

Y 

Y 

Y 

Y 

Y 

Groundwater 

transport 

1  -D  advective, 

3-D  dispersive 

1  -D  advective, 

3-D  dispersive 

1  -D  advective, 

3-D  dispersive 

1  -D  advection 
only 

1  -D  advective, 

3-D  dispersive 

1-D 

n/a  ; 

Surface-water 

transport 

1  -D  advection  & 
2-D  dispersion 

0-D  complete 
mixing 

0-D  complete 
mixing 

0-D  complete 
mixing  (onsite 
pond) 

1  -D  advection  & 
dispersion 

1-D 

Y 

Atmospheric 
transport  (far 
field) 

Gaussian  Plume 
model 

Gaussian  Plume 
mode! 

Gaussian  Plume 
model 

Gaussian  Plume 
model 

Gaussian  Plume 
model 

Gaussian  Plume 
model 

~ 

Atmospheric 
transport  (near 
field) 

N 

N 

N 

Y  (0-D  complete 
mixing) 

Y  (0-D  complete 
mixing) 

N 

— 

Unsat.  zone  flow 

1-DSSor 

Dynamic 

1-D 

1-D 

1-D 

1-D 

1-D  SS  or 
Dynamic 

N/A 

Maximum  unsat. 
layers 

5 

20 

20 

— 

— 

5 

N/A 

Sat.  zone  flow 

Constant  & 
uniform  (1  -D) 

Constant  & 
uniform  (1  -D) 

Constant  & 
uniform  (1-D) 

Constant  & 
uniform  (1  -D) 

Constant  & 
uniform  (1-D) 

Constant  & 
uniform  (1-D) 

N/A 

Maximum  sat. 
layers 

1 

1 

1 

— 

— 

1 

N/A 

Surface-water 

flow 

2-D  steady  state 

1  -D  steady  state 

1  -D  steady  state 

Surface 

hydrology 

1  -D  steady  state 

2-D  steady  state 

Steady  state 

Multi-porosity 

N 

N 

N 

N 

N 

N 

N 

Multi-pathway 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Transformation 

products 

Y 

N 

Y 

Y 

N 

Y 

Y 

Linear 

partitioning 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Non-linear 

partitioning 

N 

N 

N 

N 

N 

N 

N 

Multiple  sources 
per  simulation 

Y 

N 

N 

N 

N 

Y 

—  | 

L 

(Continued) 
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Table  3  (Concluded) 


Mode! 

Property 

MEPAS 

MULTIMED 

MULTIMED- 

DP 

RESRAD- 

CHEM 

MMSOILS 

FRAMES 

PIRANHA 

Multiple 
contaminants 
per  simulation 

Y 

N 

N 

Y 

Y 

Y 

Y 

PC  platform 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Workstation 

platform 

Y 

N 

N 

N 

— 

Y 

Y 

1  NAPL  =  Norvaqi 

jeous  phase  liquid 

-  J 

Table  4 

MEPAS  Database  Parameters  for  TNT,  RDX,  and  HMX 

Property/Parameter 

TNT  Value 

RDX  Value 

HMX  Value 

Physical  Properties 

Molecular  weight  (g/mole) 

227.1 

222.1 

296.2 

Vapor  pressure  (mm  HG  @  25  °C) 

8.99E-08 

4.08E-04 

8.40E-04 

Henry’s  Law  constant  (atm  m3/mole  @  25  °C) 

1.37E-04 

8.30E-06 

3.53E-06 

Water  solubility  (mg/L  @  25  °C) 

480.00 

38.40 

2.63 

Carbon  matter  partition  coefficient  (ml/g) 

17.80 

1.01E+03 

2.40 

Octanol-water  partition  coefficient  (ml/g) 

77.60 

1.91E+03 

3.90 

Aqueous  skin  permeability  (cm/hr) 

0.00 

0.00 

0.00 

Environmental  Half-times 

In  air  (days) 

0.00 

0.00 

0.00 

In  groundwater  (days) 

0.00 

0.00 

0.00 

In  surface  water  (days) 

5.78E-01 

6.63 

17.70 

In  soil  (days) 

0.00 

2.52E+06 

1.17E+06 

Environmental  Transfer  Factors  5 

Bioaccumulation  in  fish  (LVkg) 

0.00 

0.00 

0.00 

Bioaccumulation  in  shellfish  (L/kg) 

0.00 

0.00 

0.00 

Soil-to-plant  transfer  factor 

0.00 

0.00 

0.00 

Feed-to-animal  meat  transfer  factor  (d/kg) 

0.00 

0.00 

0.00 

Feed-to-cow  milk  transfer  factor  (d/L) 

0.00 

0.00 

0.00 

Grout  diffusion  coefficient  (cm2/sec) 

5.00E-08 

5.00E-08 

5.00E-08 

Water  purification  factor  (mg/kg/day) 

1.00 

1.00 

1.00 

Deposition  velocity  (m/sec) 

0.00 

0.00 

0.00 

Atmospheric  deposition  class 

6 

6 

6 

Toxiclties  and  Slope  Factors 

■  Jt  ' 

Inhalation  slope  factor  (mg/kg/day) 

3.00E-02 

0.11 

0.00 

Ingestion  slope  factor  (mg/kg/day) 

3.00E-02 

0.11 

0.00 

Inhalation  reference  dose  (mg/kg/day) 

5.00E-04 

3.00E-02 

5.00E-02 

Ingestion  reference  dose  (mg/kg/day) 

5.00E-04 

3.00E-02 

5.00E-02 

Gastro-intestinal  absorption  factor 

0.00 

0.00 

0.00 

Skin  absorption  fraction  from  soil 

0.00 

0.00 

0.00 
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Proposed  revisions  to  the  MEPAS  database  for  TNT,  RDX,  and 
HMX 

Evaluation  of  some  of  the  most  current  data  and  analyses  available  led 
to  recommended  changes  to  parameters  in  the  MEPAS  database.  Table  5 
shows  the  environmental  data  (Townsend  and  Myers  1996)  where  changes 
to  the  MEPAS  database  are  suggested  for  TNT,  RDX,  and  HMX. 


Table  5 

Literature  Values  for  Key  TNT,  RDX,  and  HMX  Environmental  Data 

...... 

Property 

TNT  Value 

RDX  Value 

HMX  Value 

Vapor  pressure  (mm  Hg  @  25  °C) 

5.51  E-06 

4.03E-09 

3.33E-14 

Henry’s  Law  constant  (atm  m3/mole  @  25  °C) 

1.10E-08 

1.96E-11 

2.60E-1 5 

Water  solubility  (mg/L  @  25  °C) 

100.00 

45.00 

5.00 

Carbon  matter  partition  coefficient  (ml/g) 

524.81 

7.76- 269.1 51 

3.47  -  676.081 

Octanol-water  partition  coefficient  (ml/g) 

72.4-114.81 

CM 

1.15-1.82' 

Environmental  half-time  in  groundwater  (days) 

28  - 3601 

| 1  Denotes  range  of  values  cited  in  literature  (Townsend  and  Myers  1 996)  j| 

Brannon  et  al.  (in  preparation)  studied  the  behavior  of  TNT  and  its 
mono  and  diamino  transformation  products  in  a  surface  and  aquifer  soil 
under  both  aerobic  and  anoxic  conditions  and  determined  the  sorption 
coefficients  for  these  compounds.  Table  6  shows  the  distribution  coeffi¬ 
cients  (also  called  adsorption  coefficients)  obtained  for  TNT  and  the  trans¬ 
formation  products  4ADNT,  2  AD  NT,  2,4-DANT,  and  2,6-DANT  under 
aerobic  conditions  for  a  surface  and  an  aquifer  soil.  The  properties  of  the 
soils  which  were  used  in  the  Brannon  study  are  shown  in  Table  7. 


Table  6 

Distribution  Coefficients  (Kd)  for  TNT  and  its  Transformation 
Products  Obtained  Under  Aerobic  Conditions  in  a  Surface  and  an 
Aquifer  Soil  (from  Brannon  et  al.  in  preparation) 

Compound 

Sharkey  Clay  Surface  Soil 

LAAP  Aquifer  Soil 

Kd,  ml/g 

Std.  Error 

Kd,  ml/g 

Std.  Error 

- 

TNT 

2.2 

0.12 

0.99 

0.23 

0.02 

0.99 

4  AD  NT 

4.0 

0.22 

0.99 

0.22 

0.04 

0.86 

2  AD  NT 

5.5 

0.26 

0.99 

0.31 

0.05 

0.92 

2,4-DANT 

2.1 

0.27 

0.94 

0.08 

0.02 

0.79 

2,6-DANT 

5.5 

0.48 

0.97 

0.09 

0.002 

0.99 
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Table  7 

Soil  Properties  (from  Brannon  et  al.  in  preparation) 

Soil 

TOC1,  % 

%  Sand 

%  Silt 

%  Clay 

CEC 

meq/1  OOg2 

Sharkey  clay 

2.4 

13.8 

37.5 

48.7 

38.9 

LAAP  SP-SM 
Sandy  silt 

0.015 

92.5 

2.5 

5.0 

3.6 

1  Total  organic  carbon 

2  Cation  exchange  capacity  expressed  in  milliequivalents  per  100  grams  soil 

Equation  1  shows  the  theoretical  linear,  reversible  equilibrium  parti¬ 
tioning  coefficient  relationship 


Kd  = 


Q 


(i) 


where 

3 

Kd  =  distribution  coefficient,  L  /M 

C  =  mass  of  chemical  sorbed  per  mass  of  solid,  M/M 

Cd  =  dissolved  chemical  mass  per  unit  volume  of  liquid,  M/L 

M  =  mass 

L  =  length 

The  MEPAS  model  allows  the  user  four  options  when  specifying  the 
adsorption  coefficient  Kd.  If  the  adsorption  coefficient  is  known  to  the 
user,  then  that  value  may  be  used  directly  within  MEPAS  from  the  “Work¬ 
sheet”  menu  and  by  selecting  template  2.7-“Adsorption  Coefficients 
(. Kd)”.  From  this  worksheet,  the  user  may  directly  enter  the  adsorption 
coefficient.  The  second  option  is  applicable  if  the  user  does  not  know  the 
adsorption  coefficient  but  does  know  the  fraction  of  organic  carbon  (foc ) 
and  the  organic  carbon  partition  coefficient  (Koc).  In  this  case,  Equation  2 
(Karickhoff  1979)  may  be  used  for  most  nonpolar  organics  to  determine 
the  adsorption  coefficient  (Kd). 


Kd  =  foc*Koc  (2) 

where 

foc  -  fraction  of  organic  carbon 
Kqc  -  organic  carbon  partition  coefficient,  ml/g. 
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The  third  option  applies  if  the  user  knows  the  foc  but  does  not  know  the 
Koc.  In  this  case,  the  user  can  estimate  Kd  from  Equation  2  but  use  the 
value  of  Koc  obtained  from  the  MEPAS  environmental  database  for  the 
chemical  constituent  of  concern.  The  Koc  is  readily  obtainable  for  a  par¬ 
ticular  constituent  from  the  MEPAS  environmental  database.  The  com¬ 
plete  listing  of  data  for  a  constituent  is  obtained  by  selecting  the  “Site” 
option  from  the  MEPAS  main  menu,  then  selecting  the  “Edit  Site”  option. 
This  brings  up  the  “Site  Summary”  menu  where  the  user  may  select  the 
“Constituent  Source  Defined”  option  and  press  Enter  twice.  The  fourth 
option  is  to  use  the  default  values  which  MEPAS  provides  in  the  case 
where  the  user  does  not  know  the  foc  nor  the  Koc.  In  this  case,  MEPAS 
will  calculate  a  suggested  value  based  on  the  user’s  description  of  the  soil 
properties  (percentage  of  sand,  silt,  clay  and  organic  matter)  and  the  fol¬ 
lowing  relationship  between  the  soil  properties  and  the  foc  as  proposed  by 
Whelan  et  al.  (1987): 


foc  =  (0.005  *  %  sand  +  0.4  *  %  silt 

+ 2.0  *%  clay  +  57.735  *  %  om)  1IT4 


(3) 


where 

%sand  =  percentage  of  sand  component  of  soil 
%silt  =  percentage  of  silt  component  of  soil 
%clay  =  percentage  of  clay  component  of  soil 
%om  =  percentage  of  organic  matter  in  the  soil 

In  order  to  provide  the  user  with  a  more  realistic  estimated  value  of  Kd 
for  TNT,  the  MEPAS  database  was  updated  using  the  information  in 
Tables  6  and  7  from  Brannon  et  al.  (in  preparation)  for  TNT  in  Sharkey 
clay.  Rearranging  Equation  2  to  solve  for  the  value  Koc  yields  Equation  4: 


Joc  (4) 

When  Equation  4  is  solved,  it  provides  a  Koc  value  of  91.67  ml/g.  The 
Koc  value  obtained  in  this  manner  is  a  factor  of  5.0  of  the  MEPAS  data¬ 
base  original  default  value  for  TNT  of  17.80  ml/g.  The  original  Koc  value 
in  the  MEPAS  database  was  referenced  to  the  USEPA’s  GEMS  VP-1, 2,K 
methods. 

The  data  from  Tables  6  and  7  for  LAAP  aquifer  soil  were  not  used  in 
the  calculation  of  Koc  from  Equation  4  because  the  fraction  of  organic 
carbon  in  the  aquifer  soil  is  less  than  0.5  percent.  Typically,  nonlinearity 
occurs  on  a  plot  of  foc  versus  Kd  for  foc  values  less  than  0.5  percent 
(Rebhun  et  al.  1992).  The/oc  formulation  in  MEPAS  is  essentially 
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unaffected  by  this  nonlinearity  because  it  uses  Equation  3  to  obtain  a 
value  for  foc,  which  allows  the  fraction  of  organic  matter  to  be  weighted 
among  the  other  components  of  sand,  silt,  and  clay. 

With  a  new  value  of  Koc  for  TNT,  a  revised  value  for  the  octanol-water 
partition  coefficient  (Kow)  can  now  be  derived  for  the  MEPAS  database. 
Karickhoff  (1984)  noted  that  predictions  of  Koc  based  on  values  of  Kow 
are  much  more  reliable  than  those  based  on  solubility  (McGrath  1995). 
The  Multimedia  Modeling  Environmental  Database  Editor  (MMEDE),  a 
database  support  program  supplied  with  MEPAS,  can  estimate  values  of 
Koc  based  on  solubility  or  Kow.  MMEDE  contains  Koc  estimation  routines 
for  pesticides,  aromatics,  assorted  chemicals,  and  s-triazines  for  a  total  of 
seven  methods  of  estimation  (Warren  and  Strenge  1994).  For  the  aromatic 
compound  TNT,  one  of  two  methods  can  be  chosen  to  estimate  Koc.  The 
first  is  the  GEMS  Method  K-2  (also  termed  the  Lyman  4-10  method) 
which  is  as  follows: 


log  Koc=  log  i^-0.21 


(5) 


where 


Kow  =  octanol-water  partition  coefficient,  ml/g. 
The  other  method  is  the  Lyman  Method  4-9 


log  Koc  =  0.937  *  (log  Kow)~  0.006 


(6) 


Rearranging  Equation  6  above  and  solving  for  Kow  based  on  a  Koc 
value  of  91.67  ml/g,  yields  a  value  for  Kow  of  126.05  ml/g.  This  value  is 
slightly  out  of  the  range  of  72.4  -114.8  ml/g  compiled  and  presented  by 
Townsend  and  Myers  (1996).  For  comparison,  solving  the  Lyman  4-10 
equation  for  Kow  yields  a  value  of  148.67  ml/g. 

Table  8  summarizes  the  revisions  made  to  the  MEPAS  database  for 
TNT.  There  are  no  published  recent  studies  with  estimated  or  measured 
partition  coefficients  for  RDX  and  HMX.  There  are  some  studies  with 
RDX  and  HMX  which  are  currently  on-going  and  with  preliminary  results 
of  Kd.  Therefore,  Table  8  shows  the  values  that  were  changed  in  the 
MEPAS  database  using  the  data  from  Townsend  and  Myers  (1996)  for 
RDX  and  HMX. 
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3  Process  Descriptions  for 
the  Fate  of  Target  MRCs 


Process  descriptions  for  explosive  compounds  are  poorly  developed 
because  specific  reaction  mechanisms  and  their  interrelations  are  poorly 
understood  (McGrath  1995).  This  section  describes  the  processes  affect¬ 
ing  the  fate  of  military  explosive  compounds  and  radionuclides  in  the  envi¬ 
ronment.  This  section  will  concentrate  more  on  explosive  rather  than 
radioactive  compounds,  since  the  processes  affecting  explosive  com¬ 
pounds  in  the  environment  have  not  been  well  studied  nor  analyzed  in  as 
great  detail  as  radioactive  compounds.  Processes  affecting  groundwater 
transport  of  explosives  include,  but  are  not  limited  to,  advection,  hydrody¬ 
namic  dispersion,  adsorption,  biodegradation,  abiotic  transformations, 
solubility,  volatilization,  and  chemical  reactions. 


Factors  Affecting  Contaminant  Transport 

Darcy  s  law  states  that  the  rate  of  flow  through  a  porous  medium  is 
directly  proportional  to  the  difference  in  head  and  inversely  proportional 
to  the  length  of  the  flow  path.  The  mathematical  representation  of  Darcy’s 
law  is 


where 

v  =  specific  discharge,  L/t 
k  =  intrinsic  permeability,  L2 
p  =  fluid  density,  M/L3 
g  =  acceleration  due  to  gravity,  L/t2 
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-1  -1 

(x  =  dynamic  viscosity,  ML  t 
dh/dl  =  hydraulic  gradient,  L/L 
L  =  length 
t  =  time 
M  =  mass 

If  one  desires  to  use  the  groundwater  terminology  where  hydraulic  con¬ 
ductivity  is  used  instead  of  permeability,  then  Equation  7  can  be  rewritten 
using  a  term  for  hydraulic  conductivity  as 


dh 

v -kTi 


where 


K  = 


kpg 


(8) 


(9) 


and 


K  =  hydraulic  conductivity,  Lit 

The  general  partial  differential  advection-dispersion  equation  (Freeze 
and  Cherry  1979)  describing  three-dimensional  solute  transport  in  a 
porous,  saturated,  homogeneous,  and  isotropic  medium  under  steady-state 
flow,  and  assuming  that  Darcy’s  law  applies,  is 


dC 

dt 


(10) 


where 


Dx 

D. 


D. 


v. 


3 

mass  of  solute  per  unit  volume  of  solution,  M/L 

2 

dispersion  coefficient  in  the  x  direction,  L  / 1 

2 

dispersion  coefficient  in  the  y  direction,  L  /t 
dispersion  coefficient  in  the  z  direction,  L  / 1 

average  linear  velocity  in  the  x  direction,  L/t 
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V 

y  =  average  linear  velocity  in  the  y  direction,  L/t 

vz  =  average  linear  velocity  in  the  z  direction,  L/t 
pb  =  bulk  mass  density  of  the  porous  medium,  M/L 
n  =  porosity 

S  =  the  mass  of  chemical  constituent  adsorbed  on  a  unit  mass  of 
solid  of  the  porous  medium,  M/M 

This  equation  is  the  basic  fundamental  principle  to  all  transport  model¬ 
ing  efforts.  The  source  term  is  added  to  the  equation  to  handle  the  effects 
of  retardation  of  solute  transport  through  adsorption,  chemical  reaction, 
biological  transformation,  or  radioactive  decay. 


Advection 

Advection  is  the  process  by  which  a  compound  is  transported  primarily 
due  to  the  bulk  motion  of  flowing  groundwater.  In  Equation  10,  that  por¬ 
tion  of  transport  due  to  advection  is 


Advection  is  applicable  to  the  transport  of  both  military  energetic  com¬ 
pounds  and  radioactive  materials. 


Hydrodynamic  dispersion 

Hydrodynamic  dispersion  is  the  process  by  which  the  contaminant  is 
transported  by  variations  in  the  microscopic  velocity  within  each  pore 
channel  and  from  one  channel  to  the  next.  It  is  caused  by  the  mechanical 
mixing  during  fluid  advection  and  the  molecular  diffusion  due  to  the 
thermal-kinetic  energy  of  the  solute  particles.  Dispersion  caused  by  the 
mechanical  mixing  of  a  solute  due  to  fluid  movement  is  called  mechanical 
dispersion.  Molecular  diffusion  is  usually  of  concern  only  at  low  veloci¬ 
ties.  Unlike  molecular  diffusion,  mechanical  dispersion  is  usually  of 
greater  concern  because  it  has  more  pronounced  effects  on  the  transport  of 
a  solute.  In  Equation  10  that  portion  of  transport  due  to  dispersion  is 
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Dispersion  can  be  an  important  transport  mechanism  for  both  radioac¬ 
tive  and  energetic  military  compounds. 


Adsorption 

Adsorption  is  the  process  by  which  molecules  adhere  to  a  surface  with 
which  they  come  into  contact.  Adsorption  depends  on  the  type  of  clay,  the 
initial  ratio  of  adsorbed  mole  fractions,  and  the  concentration  of  solute  in 
suspension. 


Biotransformation 

Biotransformation  is  the  process  by  which  a  substance  is  modified  into 
another  substance(s)  by  the  activities  of  living  organisms  (e.g.,  bacteria 
and  fungi).  A  more  specific  form  of  biotransformation  is  biodegradation. 
Biodegradation  is  the  metabolic  breakdown  of  materials  into  simpler  com¬ 
ponents  by  living  organisms.  These  other  components  can,  at  times,  be 
more  detrimental  to  the  environment  than  the  original  material.  Biotrans¬ 
formation  of  military  specific  radioactive  material  is  not  likely  to  occur. 


Abiotic  transformation 

Abiotic  transformation  is  the  process  by  which  a  chemical  substance  in 
the  environment  is  modified  by  nonbiological  mechanisms.  This  usually 
includes  polymerization,  oxidation-reduction  reactions,  and  hydrolysis.  If 
abiotic  transformation  of  substances  at  or  near  the  surface  is  considered, 
then  this  would  include  the  process  of  photolysis. 

Polymerization  is  the  reaction  process  by  which  the  molecules  of  a 
monomer  combine  to  form  larger  molecules  with  a  molecular  weight 
greater  than  that  of  the  original  monomer,  resulting  in  a  molecule  with 
repeated  structural  units. 

The  oxidation-reduction  reaction  is  an  oxidizing  chemical  change 
where  an  element’s  positive  valence  is  increased  (electron  loss)  and  an 
accompanied  simultaneous  element’s  valence  is  reduced  (electron  gain). 
Stated  another  way,  a  reaction  in  which  one  reactant  is  oxidized  (loses 
electrons)  and  one  reactant  is  reduced  (gains  electrons). 

Hydrolysis  is  the  reaction  between  any  substance  and  the  hydrogen  and 
hydroxyl  ions  found  in  water. 

Photolysis  is  the  process  by  which  ultraviolet  energy  from  sunlight  is 
used  to  break  the  chemical  bonds  of  a  substance  and  form  a  modified 
substance(s). 
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Solubility 


Solubility  is  the  mass  of  a  chemical  substance  that  will  dissolve  in  a 
unit  volume  of  solution  under  specified  conditions  of  temperature  and 
pressure.  Solubility  is  a  critical  parameter  in  modeling  an  explosive  con¬ 
taminant’s  environmental  fate  and  transport  (McGrath  1995).  Many  explo¬ 
sive  compounds  and  their  derivatives  are  present  as  crystalline  solids  in 
the  subsurface  where  they  can  dissolve  slowly  and  leach  into 
groundwaters. 


Volatilization 

Volatilization  is  the  process  by  which  a  chemical  substance,  either  in 
solid  or  aqueous  phase,  undergoes  a  transfer  of  mass  to  the  vapor  phase 
due  to  high  temperatures  and/or  low  pressures  weakening  the  chemical 
bonds  which  hold  the  substance  together  in  its  solid  or  aqueous  phase. 
Individual  components  of  the  substance  with  low  boiling  points  or  melting 
points  usually  separate  first,  followed  by  components  of  next  higher 
boiling/melting  point  order  until  an  equilibrium  point  is  reached  where 
little  or  no  more  components  are  removed.  For  most  common  explosives, 
volatilization  is  relatively  insignificant  due  to  the  high  boiling  point  or 
melting  point  of  the  compounds  involved. 


Energetic  Military  Compounds 


Fate  of  TNT 

The  major  factors  affecting  fate  and  transport  of  TNT  in  the  subsurface 
include  abiotic  transformations,  sorption,  diffusion,  advection,  hydrody¬ 
namic  dispersion,  and  biodegradation  (McGrath  1995).  Recent  research 
indicates  that  biotic  transformations  are  even  more  pronounced  than 
abiotic  transformations  (Townsend  and  Myers  1996).  Some  of  these  proc¬ 
esses  are  more  active  and  prominent  than  others,  and  some  may  not  even 
have  a  mentionable  role  for  a  particular  system. 

a .  Transformation  of  TNT.  McCormick,  Feeherry,  and  Levinson  (1976) 
discovered  that  TNT  transforms  under  both  aerobic  and  anaerobic 
conditions.  Depending  on  the  reducing  potential  of  the  system,  one 
to  three  of  the  nitro  groups  could  be  reduced  to  amino  groups.  The 
transformation  pathway  of  TNT  was  studied  by  Kaplan  and  Kaplan 
(1982)  who  suggested  the  routes  depicted  in  Figure  1.  Several 
studies  have  shown  higher  production  of  4-amino-2,6-dinitrotoluene 
(4A-DNT)  than  2-amino-4,6-dinitrotoluene  (Olin,  Myers,  and 
Townsend  1996).  Townsend  and  Myers  (1996)  indicated  that  the 
preferred  reduction  route  is  TNT  — *  4A-DNT  2,4-diamino-6- 
nitrotoluene  (2,4-DANT). 
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Figure  1.  TNT  transformation  pathway  diagram  (McGrath  1995)  (after  Kaplan 
and  Kaplan  1982) 


TNT  transformations  in  unsterilized  soils  are  far  more  extensive  than 
those  in  sterilized  soils,  indicating  that  the  abiotic  transformation  compo¬ 
nent  is  small  compared  to  the  biotic  transformation  component  (Townsend 
and  Myers  1996).  Studies  have  indicated  that  TNT  transformations  are  sig¬ 
nificantly  enhanced  under  anaerobic  conditions  (Townsend  and  Myers 
1996). 

b.  Sorption  of  TNT.  The  determination  of  sorption  parameters  for  TNT 
is  often  made  difficult  because  of  the  formation  of  transformation 
products  such  as  4A-DNT,  2A-DNT,  2,4-DANT,  and  2,6-DANT 
(Brannon  and  Myers  1997).  Cosolute  effects  on  the  sorption  of  TNT 
were  studied  by  Ainsworth  et  al.  (1993),  and  their  data  suggest  that 
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TNT  and  RDX  utilize  at  least  a  fraction  of  the  same  sorption  sites. 
Competition  by  these  solutes  for  sorption  sites  reduces  the  sorption 
of  individual  solutes  as  compared  with  their  single  solute  sorption 
values  (Townsend  and  Myers  1996).  In  addition  to  competition 
among  cosolutes,  competition  for  sorption  sites  among  TNT  and  its 
transformation  products  occurs  (Townsend  and  Myers  1996).  Stud¬ 
ies  have  reported  distribution  coefficients  for  TNT  in  the  range  0.0 
to  56  ml/g  (Townsend  and  Myers  1996). 


Fate  of  RDX 

The  major  factors  affecting  the  fate  and  transport  of  RDX  in  the  subsur¬ 
face  are  abiotic  transformation,  advection,  sorption,  diffusion,  hydrody¬ 
namic  dispersion,  and  biodegradation. 

a.  Transformation  of  RDX.  The  transformation  of  RDX  has  not  been 
well  documented  because  of  the  lack  of  microanalytical  procedures 
for  measuring  the  transformation  products  (Townsend  and  Myers 
1996).  A  scheme  for  the  biodegradation  of  RDX  was  proposed  by 
McCormick,  Cornell,  and  Kaplan  (1981).  In  their  proposed  scheme 
(shown  in  Figure  2),  RDX  degrades  through  a  series  of  nitro  group 
reductions  to  a  point  where  ring  fragmentation  occurs.  This  study 
found  that  RDX  degradation  occurred  only  under  anaerobic  condi¬ 
tions;  however,  a  recent  study  by  Myers  et  al.  (in  preparation) 
showed  RDX  degradation  occurring  under  both  aerobic  and  anaero¬ 
bic  conditions.  Reports  on  the  transformation  of  RDX  have  shown 
that  RDX  transformation  is  much  less  significant  than  TNT 
transformation. 

b.  Sorption  of  RDX.  In  a  study  by  Ainsworth  et  al.  (1993),  cosolute 
effects  on  the  sorption  of  RDX  suggest  that  RDX  and  TNT  compete 
for  the  same  sorption  sites.  Distribution  coefficients  for  RDX  have 
been  reported  as  being  measured  in  the  range  from  0.0  to  6.75  ml/g 
(Townsend  and  Myers  1996). 


Transformation  of  HMX 

The  transformation  of  HMX  has  little  study,  but  those  studies  that  have 
been  performed  reported  distribution  coefficients  in  the  range  0.0  to 
13.25  ml/g  (Townsend  and  Myers  1996).  Like  RDX,  HMX  transformation 
is  much  less  significant  than  that  of  TNT. 
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Figure  2.  RDX  biodegradation  scheme  (McGrath  1995)  (after  McCormick,  Cornell,  and  Kaplan  1981) 
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Radioactivity 


Radioactive  decay  process 

Radioactive  decay  is  an  exponential  rate  process  in  which  the  number 
of  atoms  likely  to  decay  in  a  given  amount  of  time  is  proportional  to  the 
number  of  atoms  present.  This  is  expressed  mathematically  as  a  first- 
order  differential  equation  in  the  form, 


dN 

dt 


kN 


(11) 


where 

N  =  number  of  atoms 

dN/dt  =  instantaneous  time  rate  of  change  of  the  number  of  atoms 
k  =  decay  constant 

Integrating  Equation  11  above  will  provide  the  number  of  atoms  at  any 
time 


N{t)=N0e-la  (12 

where 

N(t)  =  number  of  atoms  at  time  t 
Nq  =  number  of  atoms  present  at  time  equal  zero 
e  =  base  of  the  natural  logarithm 
k  =  decay  constant 
t  =  time 

Equation  12  can  also  be  expressed  in  terms  of  half-life,  which  is  the 
amount  of  time  required  for  one-half  of  the  number  of  atom  nuclei  to 
decay.  The  relation  between  the  decay  constant  and  half-life  is  as  follows, 


log,  2 


(13) 


where 


k  =  decay  constant 
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ty2  =  half-life 


Using  the  relation  between  the  decay  constant  and  the  half-life,  Equa¬ 
tion  12  can  be  rewritten  in  the  following  form 


N(t)  =  N0(05)‘i 

where 

N(t)  =  number  of  atoms  at  time  t 
Nq  =  initial  number  of  atoms  at  time  equal  zero 
t  =  time 
ty2  =  half-life 


(14) 


Radioactive  modes  of  decay 

Radioactive  decay  occurs  through  four  primary  decay  modes.  These  pri¬ 
mary  radioactive  decay  modes  include  alpha  decay,  beta  decay,  positron 
decay,  and  electron  capture.  Each  of  these  types  of  radioactive  decay  can 
create  a  newly  formed  nucleus  with  an  excess  of  energy  which  is  usually 
released  by  emitting  one  or  more  gamma  rays.  The  measurement  of  the 
gamma  rays  is  useful  in  identifying  the  type  of  radioactivity  present  since 
the  gamma  rays  have  high  energy  which  is  characteristic  of  the  particular 
radioactive  nucleus. 


a.  Alpha  decay.  Alpha  particles  are  stable  entities,  identical  to  a 
helium-4  nucleus,  which  consist  of  two  protons  and  two  neutrons. 
Alpha  decay  is  a  type  of  radioactive  disintegration  caused  when  an 
unstable  atomic  nucleus  dissipates  excess  energy  by  spontaneously 
ejecting  an  alpha  particle.  The  following  is  an  example  of  alpha 
decay: 


241 


Am 


237 


Np  +  a 


where  americium-241  decays  to  form  neptunium-237,  and  a  represents 
the  release  of  an  alpha  particle. 

b.  Beta  decay.  When  the  nucleus  of  an  atom  has  too  many  neutrons,  the 
most  likely  course  of  decay  is  that  the  nucleus  will  emit  an  electron. 
Historically,  electrons  that  are  released  in  this  manner  are  referred 
to  as  beta  particles  or  rays,  having  been  named  before  they  were 
identified  as  electrons.  The  following  is  an  example  of  beta  decay: 


239 


Pu  +  P 
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where  neptunium-239  decays  into  plutonium-239,  and  p  represents  the 
release  of  a  beta  particle  (or  electron). 

c.  Positron  decay.  When  the  nucleus  of  an  atom  has  too  few  neutrons  it 
may  acquire  another  neutron  by  ejecting  a  positive  electron,  or  posi¬ 
tron.  This  relieves  the  nucleus  of  one  unit  of  positive  charge.  An 
example  of  this  kind  of  radioactive  decay  is  the  following: 


ll. 


ll 


B  +  p+ 


where  carbon-11  decays  into  boron-11,  and  p+  represents  the  release  of  a 
positron.  The  positron  is  composed  of  antimatter  and  is  very  short-lived 
because  of  its  existence  in  a  world  composed  of  matter.  When  a  positron 
meets  an  ordinary  electron  the  two  particles  annihilate  each  other  and 
release  their  energy  as  two  gamma  rays  of  0.5  megaelectron  volts  (MeV) 
each.  Positron  radioactivity  is  always  accompanied  by  the  emission  of 
gamma  rays  of  0.5  MeV  in  addition  to  any  other  gamma  rays  which  may 
be  emitted. 

d.  Electron  capture.  Another  way  that  the  nucleus  of  an  atom  with  too 
few  neutrons  may  gain  a  neutron  is  to  capture  one  of  the  negatively 
charged  electrons  orbiting  about  the  nucleus.  An  electron  from  an 
inner  shell  of  an  atom  is  absorbed  by  the  nucleus.  In  the  process,  a 
proton  is  converted  into  a  neutron,  and  the  atomic  number  is 
decreased  by  one.  An  example  of  this  type  of  radioactive  decay  is  as 
follows: 


7  7 

7Be  +  e-  -*  Li 


where  beryllium-7  decays  to  form  lithium-7. 


Radioactive  elements  associated  with  MRCs  and  their  fate 

Some  of  the  radioactive  elements  used  or  created  in  the  manufacture  of 
nuclear  weapons  include,  but  are  not  limited  to,  such  elements  as  uranium- 
233,  234,  235,  238,  239;  neptunium-237,  239;  plutonium-238,  239,  240, 
241,  242;  americium-  241;  thorium-232;  and  tritium  (a  radioactive  isotope 
of  hydrogen). 

a.  Fate  of  plutonium.  Plutonium  is  a  silvery-white  radioactive  metal 
that  exists  as  a  solid  under  normal  conditions.  Small  amounts  of  plu¬ 
tonium  occur  naturally,  but  by-far  the  majority  has  been  produced 
by  man  in  nuclear  reactors.  Plutonium  is  produced  when  uranium 
absorbs  an  atomic  particle.  Plutonium  has  several  isotopes,  some  of 
the  most  common  being  plutonium-238  and  plutonium-239,  where 
plutonium-239  is  the  primary  source  of  nuclear-weapon-grade  pluto¬ 
nium.  Because  plutonium  is  radioactive,  it  constantly  decays, 
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releasing  energy  and  forming  new  products.  The  energy  released  is 
called  alpha  radiation  because  of  the  release  of  alpha  particles.  The 
new  products  of  the  plutonium  decay  process  are  called  daughter 
products.  Daughter  products  are  radioactive  as  well,  and  they,  too, 
continue  to  decay  until  a  nonradioactive  daughter  is  formed. 
Plutonium-238  has  a  half-life  of  approximately  90  years  and  that  of 
plutonium-239  is  24,000  years  (ATSDR  1990). 

A  review  of  the  MEPAS  database  for  data  on  plutonium  is  presented  in 
Table  9.  A  search  through  literature  for  justification  of  any  changes  to  the 
database  found  none.  The  daughter  products  for  five  of  the  plutonium  iso¬ 
topes  used  in  the  MEPAS  database  are  shown  in  Table  10. 

■2 

b.  Fate  of  tritium.  Tritium  (  H)  is  a  radioactive  isotope  of  hydrogen  and 
is  used  in  the  manufacture  of  nuclear  weapons  as  well  as  in  the 
nuclear  power  industry.  Tritium  undergoes  beta  decay  and  forms 
helium-3  (3He)  as  the  daughter  product  which  is  stable  (i.e.,  nonra¬ 
dioactive).  Tritium  is  one  of  the  weakest  beta  emitters  known,  and 
the  range  of  the  beta  particles  emitted  is  approximately  5  mm  in  air 
or  0.005  mm  in  water  and  soft  tissue  (DOE  primer  on  tritium  safe 
handling  practices).  All  chemical  reactions  involving  hydrogen  can 
also  be  performed  with  tritium,  sometimes  at  a  higher  rate  if  the  trit¬ 
ium  concentration  is  high  enough  to  catalyze  the  reaction. 

A  review  of  literature  about  tritium  provided  no  justification  for 
changes  to  the  MEPAS  database  for  tritium  environmental  parameters. 

The  MEPAS  database  information  is  shown  in  Table  11. 
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Table  9 

Existing  MEPAS  Database  Information  on  Plutonium 


Molecular  weight 
(g/mole) 


Ingestion  dose  factor 
(rem/pCi) 


Physical  Data 


Dose  Factors 


inhalation  dose  factor  2.74E-08  2.28E-12  1.45E-04  1.97E-09  3.92E-04  4  29E-04 

(rem/pCi) 

Dermal  dose  factor  1.00E-07  1.00E-09  5.60E-04  5.60E-09  3.20E-03  3.60E-03  3.60E-03  I  6  90E-05 

(rem/pCi) 


1 .97E-13 

6.03E-1 1 

1.51E-13 

1.28E-13 

Environmental  Transfer  Factors 


Bioaccumulation  in  1  .OOE+02 
shellfish  (L/kg) 


2.50E+02 


1 .00E+02  1 .00E+02  1 .00E+02  I  1 .00E+02 


Soil-to-plant  transfer  5.00E-05  5.00E-05  5.00E-05  5.00E-05  5.00E-05  5  OOE-05 

factor  (d/kg) 
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Table  10 

MEPAS  Plutonium  Decay  Products 

Isotope 

Decay 

Product  1 

Decay 

Product  2 

Decay 

Product  3 

Decay 

Product  4 

Decay 

Product  5 

238Pu 

234u 

230rh 

226Ra 

222Rn 

210Pb 

239pu 

235u 

231  Pa 

227 he 

227Th 

^Ra 

240Pu 

236u 

232Th 

241  Pu 

241  Am 

237Np 

233 Pa 

233u 

242Pu 

233  u 

234Th 

234u 

230Th 

^Ra 

Table  11 

Existing  MEPAS  Database  Information  on  Tritium 

Property 

Value 

Physical  Data 

Molecular  weight  (g/mole) 

20 

Vapor  pressure  (mm  Hg) 

2.38E+01 

Henry’s  Law  constant  (atm  m3/mole) 

6.00E-07 

Half-life  (days) 

4.49E+03 

Dose  Factors 

Ingestion  dose  factor  (rem/pCi) 

6.30E-1 1 

Inhalation  dose  factor  (rem/pCi) 

6.30E-11 

Dermal  dose  factor  (rem/pCi) 

6.30E-11 

Air  immersion  factor  (rem/hr  per  pCi/m3) 

0.0 

Ground  exposure  factor  (rem/hr  per  pCi/m2) 

0.0 

Water  immersion  factor  (rem/hr  per  pCi/L) 

0.0 

Environmental  Transfer  Functions 

Bioaccumulation  in  fish  (L/kg) 

1.00E+00 

Bioaccumulation  in  shellfish  (L/kg) 

1.00E+00 

Soil-to-plant  transfer  factor  (d/kg) 

0.0 

Feed-to-animal  meat  transfer  factor  (d/kg) 

0.0 

Feed-to-cow  milk  transfer  factor  (d/L) 

0.0 

Grout  diffusion  coefficient  (cm2/sec) 

5.00E-08 

Water  purificiation  factor 

1 .00E+00 

Deposition  velocity  (m/sec) 

0.0 

Atmospheric  deposition  class 

4 

Toxicities  and  Slope  Factors 

Inhalation  slope  factor  (risk/pCi) 

9.59E-14 

Ingestion  slope  factor  (risk/pCi) 

7.15E-14 

External  exposure  (risk/yr  per  pCi/g) 

0.0 

Gastro -intestinal  absorption  factor 

1 .00E+00 

Skin  absorption  fraction  from  soil 

0.0 
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c.  Fate  of  uranium.  Uranium  is  a  heavy,  silvery-white  metal  that  has 
16  isotopes,  all  of  which  are  radioactive.  Naturally  occurring 
uranium  nominally  contains  99.28305  percent  by  weight  238U, 
0.7110  percent  235U,  and  0.0054  percent  234U  (Hammond  1995). 
Uranium-238  is  of  particular  importance  because  it  can  be  readily 
converted  through  a  series  of  decay  reactions  into  plutonium.  Acids 
dissolve  the  metal,  but  it  is  unaffected  by  alkalis. 

A  review  of  the  MEPAS  database  for  data  on  uranium  is  shown  in 
Table  12.  The  radioactive  decay  products  of  uranium  from  the  MEPAS 
decay  database  are  shown  in  Table  13.  A  literature  search  found  no  addi¬ 
tional  data  to  warrant  any  changes  to  these  databases  for  uranium. 

d.  Fate  of  neptunium.  Neptunium  is  an  actinide  series  metal  with  sil¬ 
very  appearance  and  15  isotopes.  The  radio-isotope  237Np  is 
obtained  as  a  by-product  from  nuclear  reactors  in  the  production  of 
plutonium. 

The  MEPAS  database  contains  information  on  3  of  the  15  isotopes  of 
neptunium  which  is  presented  in  Table  14.  The  decay  products  from  the 
MEPAS  database  for  2^Np  are  shown  in  Table  15.  A  review  of  literature 
on  neptunium  found  no  justification  for  changes  to  the  MEPAS  database 
for  this  element. 
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Table  12 

Existing  MEPAS  Database  Information  on  Uranium 


Molecular  weight  (g/mole) _ 

Half-life  (days)  _ 


Ingestion  dose  factor  (rem/pCi)  _ 


Inhalation  dose  factor  (rem/pCi) 


Dermal  dose  factor  (rem/pCi) 


Air  immersion  factor  (rem/hr  per  pCi/m3) 


Ground  exposure  factor  (rem/hr  per  pCi/m2) 


Water  immersion  factor  (rem/hr  per  pCi/L) 


Physical  Data 


Dose  Factors 


1.31E-06 


1  6.59E-04 


3.00E-05 


|  1.88E-13 


1.35E-14 


4.29E-13 


|  2.89E-07 


1.35E-04 


6.20E-06 


2.17E-13 


I  2.83E-07 


1.32E-04 


6.10E-06 


1.02E-13 


2.33E-13 


235 

2.57E+1 1 


2.66E-07 


1.23E-04 


5.60E-06 


9.59E-11 


1.97E-12 


2.1  IE-10 


Environmental  Transfer  Factors 


5.00E+01 


5.00E+01 


Bioaccumulation  in  shellfish  (L/kg) 

6.00E+01 

6.00E+01 

6.00E+01 

6.00E+01 

Soil -to -plant  transfer  factor  (d/kg) 

3.50E-03 

3.50E-03 

3.50E-03 

3.50E-03 

Feed-to-animal  meat  transfer  factor  (d/kg) 


Feed-to-cow  milk  transfer  factor  (d/L) 


Grout  diffusion  coefficient  (cm2/sec) 


Water  purification  factor 


Deposition  velocity  (m/sec) 


Atmospheric  deposition  class 


2.00E-04 


6.00E-04 


2.00E-04 


6.00E-04 


2.00E-04 


6.00E-04 


2.00E-04 


6.00E-04 


236 

8.55E+09 


2.69E-07 


1.25E-04 


5.80E-06 


6.66E-14 


8.65E-15 


1.54E-13 


5.00E+01 


6.00E+01 


3.50E-03 


2.00E-04 


238 

1.63E+12 


239 

1.63E-02 


2.55E-07 

7.73E-1 1 

IHhlJ.IMUilJMI 

5.50E-06 

0.0 

4.53E-14 

2.90E-11 

7.34E-15 

6.84E-13 

1.06E-13 

6.49E-11 

|  5.00E+01 


6.00E+01 


3.50E-03 


2.00E-04 


|  5.00E+01 


6.00E+01 


1 .00E-03 
1 


1.00E-03  1.00E-03  1.' 

1  1  1 


Toxicities  and  Slope  Factors 


Inhalation  slope  factor  (risk/pCi) 

5.29E-08 

Ingestion  slope  factor  (risk/pCi) 

8.12E-1 1 

External  exposure  (risk/yr  per  pCi/g) 

3.42E-1 1 

Gastro-intestinal  absorption  factor 

0.0 

Skin  absorption  fraction  from  soil 

0.0 

1.41E-08 
4.48E-1 1 


1 .40E-08 
4.44E-1 1 


1 .06E-09 

2.14E-1 1 

5.00E-02 

5.00E-02 

1 .30E-08 
4.52E-11 


2.63E-07 


1.32E-08 
4.21  E-11 
1.72E-1 1 


1 .24E-08 
4.27E-1 1 
1.50E-11 


OO _ 

0.0 

1 .32E-07 


Table  13 

MEPAS  Uranium  Decay  Products 


Isotope 


Decay  Product  1 


229Th _ 

230Th 


Decay  Product  2  Decay  Product  3  Decay  Product  4  Decay  Product  5 
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Table  14  II 

Existing  ME  PAS  Database  Information  on  Neptunium 

Isotope  j| 

1  Property 

237Np 

238Np 

239Np  | 

Physical  Data 

-  -  V*,,  ■ 

1  Molecular  weight  (g/mole) 

237 

238 

239 

1  Half-life  (days) 

7.82E+08 

2.12E+00 

2.36E+00 

Dose  Factors  .  J| 

Ingestion  dose  factor  (rem/pCi) 

4.44E-06 

4.00E-09 

3.26E-09 

[  Inhalation  dose  factor  (rem/pCi) 

5.40E-04 

3.70E-08 

2.51  E-09 

Dermal  dose  factor  (rem/pCi) 

4.50E-03 

0.0 

4.30E-09  1 

|  Air  immersion  factor  (rem/hr  per  pCi/m3) 

1.37E-11 

3.62E-10 

1.02E-10 

|  Ground  exposure  factor  (rem/hr  per  pCi/m2) 

3.82E-13 

7.05E-12 

2.17E-12 

|  Water  immersion  factor  (rem/hr  per  pCi/L) 

3.08E-1 1 

7.83E-10 

2.27E-10 

1  Environmental  Transfer  Factors 

Bioaccumulation  in  fish  (LTkg) 

2.50E+02 

2.50E+02 

2.50E+02 

|  Bioaccumulation  in  shellfish  (L/kq) 

4.00E+02 

4.00E+02 

4.00E+02 

II  Soil-to-plant  transfer  factor  (d/kg) 

2.35E-03 

2.35E-03 

2.35E-03  j 

|  Feed-to-animal  meat  transfer  factor  (d/kq) 

5.50E-05 

2.00E-04 

2.00E-04 

|  Feed-to-cow  milk  transfer  factor  (d/L) 

5.00E-06 

6.00E-04 

6.00E-04 

Grout  diffusion  coefficient  (cm2/sec) 

5.00E-13 

5.00E-13 

5.00E-13 

Water  purification  factor 

7.00E-01 

7.00E-01 

7.00E-01 

II  Deposition  velocity  (m/sec) 

1.00E-03 

1.00E-03 

1.00E-03 

Atmospheric  deposition  class 

1 

3 

3 

II  Toxicities  and  Slope  Factors 

||  Inhalation  slope  factor  (risk/pCi) 

3.45E-08 

4.68E-12 

2.41E-12 

||  Ingestion  slope  factor  (risk/pCi) 

2.95E-10 

4.56E-12 

4.27E-12 

|  External  exposure  (risk/yr  per  pCi/q) 

7.56E-07 

1.95E-06 

2.42E-07 

Gastro-intestinal  absorption  factor 

1.00E-03 

1.00E-03 

1 -00E-03 

I  Skin  absorption  fraction  from  soil 

0.0 

0.0 

0.0 

Table  15 

j  MEPAS  Neptunii 

im  Decay  Products 

Isotope 

Decay  Product  1 

Decay  Product  2 

Decay  Product  3 

I237Np 

233pa 

233u 

229 Th 

e.  Fate  of  americium.  Americium  is  a  silvery-white  metal  of  the  acti¬ 
nide  series.  It  is  a  synthetic  chemical  element,  undetected  in  nature. 
All  isotopes^of  americium  are  unstable,  with  the  most  important  iso¬ 
tope  being  Am  which  is  readily  available  in  large  quantities  from 
the  production  of  plutonium.  Americium  reacts  with  oxygen  and 
with  hydrogen  forming  americium  dioxide  (Am02)  and  americium 
hydride  (AmH2),  respectively. 
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The  MEPAS  database  contains  data  on  eight  of  the  isotopes  of 
americium  which  are  presented  in  Table  16.  The  decay  products 
from  the  MEPAS  database  for  three  of  the  radio-isotopes  of  ameri¬ 
cium  are  shown  in  Table  17.  A  review  of  literature  on  americium 
found  no  justification  for  changes  to  the  MEPAS  database  for  this 
element. 

/.  Fate  of  thorium.  Thorium  is  silvery-white  metal  of  the  actinide  series 
with  25  known  isotopes,  all  of  which  are  unstable.  Thorium  is 
slowly  attacked  by  water  and  does  not  dissolve  readily  in  most  com¬ 
mon  acids. 

The  MEPAS  database  contains  data  on  7  of  the  25  isotopes  of 
thorium,  presented  in  Table  18.  The  decay  products  from  the 
MEPAS  database  for  six  of  the  radio-isotopes  of  thorium  are  shown 
in  Table  19.  A  review  of  literature  on  thorium  found  no  justification 
for  changes  to  the  MEPAS  database  for  this  element. 
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Table  16 

Existing  MEPAS  Database 


Information  on  Americium 


Physical  Data 


Molecular  weight  (g/mole) 


Half-life  (days) 


Ingestion  dose  factor  (rem/pCi) 

Inhalation  dose  factor  (rem/pCi) 

Dermal  dose  factor  (rem/pCi) 

Air  immersion  factor  (rem/hr  per  pCi/m3) 
Ground  exposure  factor  (rem/hr  per  pCi/m2) 
Water  immersion  factor  (rem/hr  per  pCi/L) 


Bioaccumulation  in  fish  (LVkg) 


Bioaccumulation  in  shellfish  (LVkg) 


Soil-to-plant  transfer  factor  (d/kg) 
Feed-to-animal  meat  transfer  factor  (d/kg) 


Feed-to-cow  milk  transfer  factor  (d/L) 


Grout  diffusion  coefficient  (cm2/sec) 
Water  purification  factor 
Deposition  velocity  (m/sec) 
Atmospheric  deposition  class 

Inhalation  slope  factor  (risk/pCi) 
Ingestion  slope  factor  (risk/pCi) 


External  exposure  (risk/yr  per  pCi/g) 


Gastro-intestinal  absorption  factor 


Skin  absorption  fraction  from  soil 


5.07E-02  6.81  E-02 


Dose  Factors 

6.59E-1 1  1.32E-10  9.88E-10  2.53E-09 

2.39E-11  8.58E-10  4.59E-10  1 .84E-09 

1.00E-08  4.00E-08  1 .00E-07  8.00E-07 

2.27E-10  5.77E-10  1.38E-10  6.66E-10 

4.72E-12  1.13E-11  2.46E-12  1.31E-11 

4.98E-10  1 .25E-09  3.07E-10  1.44E-09 

Environmental  Transfer  Factors 


2.50E+02  2.50E+02  2.50E+02  2.50E+02 


1.00E+03  1.00E+03  1.00E+03  1  .QQE+03 

1.03E-04  1.03E-04  1.03E-04  1.03E-Q4 

3.50E-06  3.50E-06  3.50E-06  3.50E-06 

4.00E-07  4.00E-07  4.00E-07  4.00E-07 


1 .00E-03  1 .00E-03  1  .OOE-03  1 .00E-03 


Toxicities  and  Slope  Factors 


241 


1.58E+05 


6.70E-01  5.55E+04  2.70E+06 


2.53E-09 

3.64E-06 

1.41  E-09 

1 .84E-09 

4.44E-04 

5.85E-08 

8.00E-07 

3.70E-03 

3.30E-07 

6.66E-10 

1.09E-11 

8.19E-12 

1.31E-11 

3.65E-13 

2.08E-13 

1.44E-09 

2.50E-1 1 

1.83E-11 

2.50E+02 

2.50E+02 

2.50E+02 

2.50E+02 

1.00E+03 

1.00E+03 

1.00E+03 

^09 

1.03E-04 

1 .03E-04 

1.03E-04 

3.50E-06 

3.50E-06 

3.50E-06 

3.50E-06 

4.00E-07 

4.00E-07 

4.00E-07 

4.00E-07 

5.00E-13 

5.00E-13 

5.00E-13 

5.00E-13 

7.00E-01 

7.00E-01 

7.00E-01 

7.00E-01 

1. 00E-03 

1. OOE-03 

1. 00E-03 

1. 00E-03 

1 

1 

1 

1 

3.85E-08 

1.04E-11 

3.49E-08 

3.28E-10 

1.47E-12 

2.92E-1 0 

4.59E-09 

5.76E-09 

9. 45  E-09 

1. 00E-03 

1 .00E-03 

1. 00E-03 

0.0 

0.0 

0.0 

Table  17 

MEPAS  Americium  Decay  Products 


Isotope 


Decay  Product  1  Decay  Product  2  Decay  Product  3  Decay  Product  4  Decay  Product  5 
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Table  18 

Existing  MEPAS  Database  Information  on  Thorium 

Isotope 

Property 

227Th 

“"Th 

229Th 

230Th 

231  Th 

232Th 

234Th 

Physical  Data 

Molecular  weight  (g/mole) 

227 

228 

229 

230 

231 

232 

234 

Half-life  (days) 

1 .87E+01 

6.99E+02 

2.68E+06 

2.81  E+07 

1.06E+00 

1.41E-10 

2.41  E+01 

Dose  Factors 

Ingestion  dose  factor  (rem/pCi) 

3.81  E-08 

3.96E-07 

3.53E-06 

5.48E-07 

1.35E-09 

2.73E-06 

1.37E-08 

Inhalation  dose  factor  (rem/pCi) 

1.62E-05 

3. 42  E -04 

2.15E-03 

3.26E-04 

8.77E-10 

1.64E-03 

3.50E-08 

Dermal  dose  factor  (rem/pCi) 

5.70E-05 

1.80E-03 

1.70E-02 

2.50E-03 

6.80E-10 

1.40E-02 

3.90E-08 

Air  immersion  factor  (rem/hr  per  pCi/m3) 

6.51  E-t  1 

1.23E-12 

5.09E-1 1 

2.33E-13 

6.94E-12 

1.16E-13 

1.65E-11 

Ground  exposure  factor  (rem/hr  per  pCi/m2) 

1.38E-12 

3.13E-14 

1.14E-12 

9.99E-15 

2.47E-13 

7.34E-15 

3.64E-13 

Water  immersion  factor  (rem/hr  per  pCi/L) 

1.43E-10 

2.73E-12 

1.14E-10 

5.24E-1 3 

1.57E-11 

2.64E-13 

3.58E-11 

Environmental  Transfer  Factors 

Bioaccumulation  in  fish  (17kg) 

1.00E+02 

1.00E+02 

1.00E+02 

1.00E+02 

1.00E+02 

1 .00E+02 

1 .00E+02  | 

Bioaccumulation  in  shellfish  (17kg) 

5.00E+02 

5.00E+02 

5.00E+02 

5.00E+02 

5.00E+02 

5.00E+02 

5.00E+02 

Soi!-to-plant  transfer  factor  (d/kg) 

3.00E-05 

3.00E-05 

3.00E-05 

3.00E-05 

3.00E-05 

Feed-to-animal  meat  transfer  factor  (d/kg) 

6.00E-06 

6.00E-06 

6.00E-06 

6.00E-06 

6.00E-06 

6.00E-06 

6.00E-06 

Feed-to-cow  milk  transfer  factor  (d/L) 

5.00E-06 

5.00E-06 

5.00E-06 

5.00E-06 

5.00E-06 

5.00E-06 

5.00E-06 

Grout  diffusion  coefficient  (cm2/sec) 

1.00E-12 

1.00E-12 

1.00E-12 

1.00E-12 

0.0 

1.00E-12 

1.00E-12 

Water  purification  factor 

7.00E-01 

7.00E-01 

7.00E-01 

7.00E-01  | 

0.0 

7.00E-01 

7.00E-01 

Deposition  velocity  (m/sec) 

1.00E-03 

1 .00E-03 

1  .OOE-03 

1. 00E-03 

Atmospheric  deposition  class 

1 

1 

1 

1 

(■■■I 

{■■HI 

[■HI 

Toxicities  and  Slope  Factors 

Inhalation  slope  factor  (risk/pCi) 

4.31  E-09 

9.45E-08 

7.60E-08 

1 .72  E-08 

1 .93  E-08 

1.90E-11 

Ingestion  slope  factor  (risk/pCi) 

4.04E-1 1 

6.29E-1 1 

5.65E-1 1 

3.75E-1 1 

1.79E-12 

3.28E-11 

1.93E-11 

External  exposure  (risk/yr  per  pCi/g) 

1 .74E-07 

5.28E-06 

5.94E-08 

4.40E-1 1 

2.09E-09 

1.97E-11 

5.71  E-08 

Gastro-intestinal  absorption  factor 

2.00E-04 

2.00E-04 

2.00E-04 

2.00E-04  | 

0.0 

2.00E-04 

2.00E-04 

Skin  absorption  fraction  from  soil 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

Table  19 

MEPAS  Thorium  Decay  Products 

Isotope 

Decay  Product  1 

Decay  Product  2 

Decay  Product  3  Decay  Product  4  Decay  Product  5 

^Th 

^Ra 

^h 

224  Ra 

225Ra 

225Ac 

230Th 

^Ra 

210Pb 

210Po 

232Th 

^Ra 

^Ra 

234Th 

234u 

23ftTh 

226Ra 

222Rn 

210Pb 
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4  Summary  and 
Recommendations 


A  review  of  existing  exposure  assessment  models  was  performed,  and  a 
multimedia  exposure  assessment  model  for  evaluating  ecological  risk  of 
exposure  to  military-related  compounds  at  military  sites  has  been  devel¬ 
oped.  The  MEPAS  model  developed  by  Pacific  Northwest  Laboratory  for 
the  U.S.  DOE  was  selected,  with  changes  made  to  its  environmental 
database.  The  model  is  applicable  to  TNT,  HMX,  RDX,  and  several  radio¬ 
active  contaminants  such  as  plutonium  and  americium.  The  model  is 
primarily  intended  as  a  screening-level  tool  to  aid  in  identifying  an  eco¬ 
logical  risk  from  an  MRC.  The  model  can  be  applied  to  inactive  as  well  as 
active  air  and  water  release  sites. 

The  model  itself  uses  a  simplified  contaminant-transport  and  exposure 
algorithm  encapsulated  into  a  user-friendly  shell.  The  model  is  an  inte¬ 
grated  system  of  analytical,  semi-analytical,  and  empirically  based  mathe¬ 
matical  models  that  simulate  the  release  of  contaminants  from  sources, 
migration  and  fate  of  these  contaminants  through  various  environmental 
media,  concentrations  at  designated  receptor  locations,  and  health  effects 
to  specified  sensitive  receptors  caused  by  exposure  to  these  contaminants. 

The  model  focuses  on  addressing  long-term  average  environmental  con¬ 
ditions  resulting  from  release  to  the  environment  from  point  and  area 
sources.  Inherent  in  the  design  of  the  model  is  the  lack  of  the  capability  to 
address  acute  exposures  to  environmental  contamination. 

Recommendations  for  model  improvement  include  future  studies  to 
address  toxicities  (such  as  reference  doses),  slope  factors  (such  as  inges¬ 
tion  slope  factor),  and  environmental  transfer  factors  (such  as  bioaccumu¬ 
lation).  These  areas  need  to  be  addressed  to  determine  values  that  are 
more  specific  to  ecological  units.  Most  of  the  parameters  from  these  areas 
were  not  modified  for  this  study  and  could  potentially  be  of  significant 
importance  in  determining  the  risk  of  exposure  to  an  ecological  unit  from 
MRCs.  Additionally,  future  applications  of  the  model  to  specific  ecologi¬ 
cal  units  are  needed  to  gather  a  more  complete  understanding  of  the 
model’s  strengths  and  limitations. 
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